New self-dual solutions of SU(2) Yang-Mills theory in Euclidean Schwarzschild 

space 



Ricardo A. MosncQ 

Instituto de Matemdtica, Estatistica e Computagao Cienttfica, 
Universidade Estadual de Campinas, 13083-859, Campinas, SP, Brazil. 

Gustavo Marques Tavare^j] 
Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, 13083-970, Campinas SP, Brazil. 

(Dated: November 3, 2009) 

We present a systematic study of spherically symmetric self-dual solutions of SU(2) Yang- 
Mills theory on Euclidean Schwarzschild space. All the previously known solutions are recov- 
ered and a new one-parameter family of instantons is obtained. The newly found solutions 
have continuous actions and interpolate between the classic Charap and Duff instantons. 
We examine the physical properties of this family and show that it consists of dyons of unit 
(magnetic and electric) charge. 
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I. INTRODUCTION 

The study of SU(2) instantons, i.e., finite action smooth (anti) self-dual solutions of the SU(2) 
Yang-Mills (YM) equation, has been extensively carried out in Euclidean flat space, and numerous 
interesting properties of them are by now well established These solutions correspond to global 
minima of the action in each topological sector and, therefore, are associated with the leading 
terms in semiclassical approximations to the path integral. In particular, instantons are associated 
with tunneling between degenerate vacua in Yang-Mills theory at zero temperature. 

At the classical level, it is well known that these solutions are classified by ^(S* 3 ) = Z and 
that their action, which coincides with the first Pontryagin class, is restricted to integer values 0]. 
These results are related to the topology of flat space and are not expected to hold in manifolds 
with other topologies. 

One important feature of self-dual solutions in Euclidean backgrounds (not only in flat space) 
is that the associated energy-momentum tensor vanishes and thus such field configurations do not 
disturb the geometry of spacetime 0]. Hence, given a solution of Einstein's equation in vacuum 
and a self-dual solution of YM theory in such background, we immediately have a full solution 
of Einstein- Yang-Mills theory. Therefore, solutions of Einstein's equation in vacuum provide the 
most interesting manifolds (from a physical viewpoint) to study self-dual solutions of YM theory. 

Furthermore, although gravity is expected to have a weak effect on the perturbative sector of 
field theories, it might drastically change the nonperturbative sector. This is not a surprise since 
gravity can change the topology of spacetime and therefore strongly affect the configuration space 
for the fields. For example, a number of works have considered self-dual solutions of Yang-Mills 
theory in various Euclidean backgrounds and found solutions with no counterpart in flat space 
0, 0, Si Sj M, fl El • Even the usual tunneling interpretation does not hold in some spacetime 
geometries [K| . Another remarkable difference is that the action can take continuous values^as has 
been exhibited for the AdS, Schwarzschild and other spherically symmetric spacetimes 
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In this work we approach the question of determining, in a systematic way, the possible spher- 
ically symmetric smooth self-dual solutions of Yang-Mills theory in the Euclidean Schwarzschild 
background. As a result, we recover all the previously known solutions and, moreover, disclose a 
whole new one-parameter family of instantons with continuous actions. This family interpolates 
between the classic Charap and Duff's solutions 0, 0] with actions 1 and 2. We examine the 
physical properties of this family and show that it consists of dyons of unit (magnetic and elec- 
tric) charge. Like the previously obtained instantons in this background, these new solutions are 
associated with a constant gauge-invariant static potential (no tunneling). 

The structure of this paper is as follows. We first discuss the spherically symmetric ansatz 
employed by us, paying special attention to regularity issues that naturally arise in this context. We 
then prove the existence of the aforementioned solutions, study their global behavior and examine 
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their physical properties. Finally, we employ the recently introduced mapping method 
obtain a new variable according to which these solutions can be analytically expressed as a power 
series with apparently infinite radius of convergence. We close with some general remarks which 
include a brief discussion on the possible mathematical relevance of our results to the study of 
instantons in asymptotically locally flat (ALF) geometries 13|. 



II. ANSATZ 



The Euclidean Schwarzschild space is a noncompact complete Riemannian manifold topolog- 



ically given by M = S 



S 2 x 



<\{0}) C M: 



and endowed with the following metric, defined on the open set 



ds 2 



^W+fl-— V dr 2 + r 2 dn 2 . 
r ) \ r J 



(1) 



In the above expression m > 0, dQ 2 is the usual round metric in S 2 , and (r, r) are polar coordinates 
on IR 2 \{0}, which take values on the intervals (2m, oo) and [0, 8irm), respectively. It is well known 
that the metric ([T]) can be uniquely extended to the whole M, and that such extension is a 
nonsingular solution of the vacuum Einstein equations, with nontrivial Euler characteristic x = 2. 

Using coordinates x^ = (r, r sin 8 cos eft, r sin 6 sin 0, r cos 6), with (#, ip) the usual angular coordi- 
nates in S 2 , we will adopt the general spherically symmetric ansatz for the gauge potential 14]: 



A = — 

2i 



X^ 1 X^ X^ X f x^ "\ ( x^ 

(t, r)—dr + P{t, r)^-dx c + (7(7-, r) - l)e abc — d — J + o(r, r)d — 



(2) 



where a a are the Pauli matrices and e a b c is the totally antisymmetric tensor, with £123 = +1. 

Because of the angular character of the time coordinate, care must be taken with the boundary 
conditions to ensure regularity of the gauge potential. Clearly all of the ansatz functions must 
be periodic in time. Also, as r — > 2m, (f> must vanish and all the other functions must become 
time-independent. 

Under a gauge transformation A = U~ 1 AU + U~ 1 dU, with U = exp(— i f(r, r)x a a a /2r), the 



form of the ansatz is unchanged (this corresponds to a [/(l)-symmetry of the problem [id . Il4( |) 
and this can be used to choose a gauge in which (3 = 0. If we are interested in time-independent 
self-dual solutions 15[ we can also set 7 = 0. The (singular) gauge transformation given by 
U = exp(— iipa^/2) exp(— i0a2/2) then leads to the ansatz [6| 



A = air) dO — + a(r) sin# dip — + (cosO dip + (p(r) dr) — -. 

2i 2% 2i 



(3) 
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Note that given a singular potential in the form of Eq. ([3]) one can turn it back to the form of 
Eq. ([2]) and get a nonsingular potential as long as we ensure that the regularity constraints are 
obeyed. 

For a potential given by the above expression the self-duality equations are: 



2?n x 



a'(r) = - 1-— a(r)<f>(r), (4a) 



r 

<t>\r) = (4b) 

It immediately follows that <fi is governed by the ordinary differential equation 

r , . d 2 d> ,dd> ndd> 

2 ( r " 2m )^l + ( r " 2m )^ " ^ + r = °> ( 5 ) 
and that a is given, in terms of (f>, by 

a 2 (r) = 1 - r 2 cj)'(r). (6) 

Notice that this imposes a restriction on the solutions of Eq. ([5]), since a needs to be a real function 
(a 2 (r) > 0) in order that A be su(2)-valued 

The Lagrangian associated with these solutions is easily calculated, 

At tr(F A *F) = [(1 - a 2 )(f>] dr A dr A dH, (7) 

8tt a 8tt z or 

and using this expression one finds that the action is given by 



If 00 
S = / tv(F A *F) = 4m [(1 - a 2 )0] . (8) 

71 " ./M 2m 



The solutions we will be interested in satisfy a(2m) = and a{r) — > as r goes to infinity (see 
section ITTTj) . In this case the action depends only on the asymptotic values of <j>(r): 

S = 4m [<j>(oo) - 4>(2m)} . (9) 

We note that this expression is not gauge invariant: its simple form is a direct consequence of the 
gauge choice associated with ansatz (|3j). 



Following the approach described in [18| (see also [10|, [19|]) one can define an Abelian field 
strength given by 

This can be identified with an electromagnetic field and, in particular we find that 

fi = ^/ (r) ^ ) (10a) 



r 



x 



*r = - 3 . (iob) 

Therefore it is clear that all configurations given by ansatz ([3]) have unit magnetic charge. The 
electric charge depends on the asymptotic behavior of (f>(r). 
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III. SOLUTIONS 
A. General discussion 

Three solutions immediately arise from Eqs. ((4|), which are the classic Charap and Duff solu- 
tions jjj: 

i. 4>{r) = 0, a(r) = 1. In this case, 5 = 0. 

ii. (j)(r) = a 2 (r) = 1 - ^. In this case, 5=1. 

iii. 4>{r) = c — -, where c is a constant, and a(r) = 0. In this case, 5 = 2. 

We see from Eqs. (jlOJl that the solutions with actions 1 and 2 above correspond to a monopole and 
a dyon, respectively [10(. The third solution is moreover Abelian, since it lies in a fixed u(l) inside 
5u(2) as can be seen from Eq. ([3]) (see also j2o|). 

It is evident that the second and third (for general c) solutions above do not satisfy the regularity 
condition <f)(2m) = 0. Going back to Eq. ([2]) we can see that given a gauge potential with <p(2m) = 
A / we can turn it into a regular potential by performing a singular gauge transformation with 
U = exp(i\t x a a a /2r), as long as 8irmA = 2irp, with p £ Z, and a(2m) = 0. This is clearly the 
case for the aforementioned solutions if we choose c properly. Notice that, unless a(r) = 0, the 
nonsingular solution will be time-dependent. 

A straightforward application of the power series method and the aforementioned condition on 
(j)(2m) = restrict the possible solutions of ©, in a neighborhood of 2m, to 4>{r) = YlT=o flfc ( r ~~ 
2m) k , with ao = -^,P € Z. It turns out that only the Abelian solution is possible for p > 0. For 
p < a much more interesting situation arises. In this case, the coefficient a„ p+ i is not fixed by the 
condition a o = and the remaining coefficients a& are determined by a complicated nonlinear 
recurrence relation that depends on all the coefficients with index less than k (and hence on a_ p+ i). 
As a result, for each p < 0, the free parameter a_ p +i gives rise to a whole family of solutions of ([5]). 

We have strong numerical evidence that, for p < — 2, the only possible solution with finite action 
is again the Abelian solution: all other possible solutions of Eq. ([5]) either diverge as r — > 00 or 
render a imaginary. In this way, we turn our attention to the two remaining possibilities, namely 
p = —1 and p = 0. These two cases have solutions with very dissimilar behaviors. In fact, as 
r — > 2m, a(r) goes like (r — 2m) 1 / 2 and const + (r — 2m) for p = —1 and p = 0, respectively. More 
to the point, the solutions associated with p = — 1 and p = are not gauge equivalent (except 
for the trivial Abelian case). This can be easily seen by plugging the power series expressions for 
these solutions into the gauge-invariant expression J7]). It is not difficult to show that the case 
p = exactly corresponds to the results obtained in [gj, in which the authors present a family of 
solutions with action ranging from to 2. Therefore, in this work we focus on the remaining possible 
choice (p = —1) which, as discussed above, leads to time-dependent solutions in the appropriate 
nonsingular gauge. 

It follows from the series solution that a\ = cf)'(2m) = j^i- Writing s = r — 2m, the first few 
terms of the corresponding solution are given by 

(f> K (s + 2m) = -— + — - + s 2 j s 3 = s 4 + 0(s 5 ), (11) 

4m 4m^ lGm^ 16m 4 256m b 

where we introduced the free dimensionless parameter k given by k = 16m 3 02- It is not hard to 
show, directly from the recurrence relation satisfied by a n , that: 

1. for k = —3, <j) K reduces to </>_3(r) = — which is the aforementioned Charap and Duff's 
solution with 5=1; 
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2. for re = —2, (j) K reduces to </>_2(r) 
solution with S = 2. 



l 

4m 



which is the aforementioned Charap and Duff's 



Unfortunately, the radius of convergence of (jlip is, in general, not greater than 2m. In this way, 
this method fails to provide any insight about (j){r) for r>l. In particular, it yields no information 
whatsoever on the finiteness (not to say the value) of the action associated with each <p K (Eq. ([9])). 
This kind of question requires considerations on the global behavior of these solutions and will be 
elucidated in the last part of this section by studying some analytic properties shared by them. 
Then, in section IV, we will use a mapping method to obtain a series solution that has an infinite 
radius of convergence. 

Despite its aforementioned limitations, Eq. can be used to generate initial values 

(4> K (2m + e), 4>' K (2m + e)) for ([5]), with e positive and arbitrarily small. Such initial values can then 
be employed to numerically integrate ([5]) [21| . Figure [T] shows representative solutions obtained in 
this way. The thick curves correspond to the classic Charap and Duff solutions Q|. 



■Q- 




r/m 

FIG. 1: Plots of 4> K {r) for evenly spaced values of re going from —4 (lowest curve) to —1 (uppermost curve). 
The thick curves correspond to the Charap and Duff solutions (re = —3 for the lower thick curve and re = —2 
for the upper thick curve). 



Figure Q] suggests that: 

1. solutions with re < —3 (below the lower thick curve) are not limited, and therefore do not 
have finite action (cf Eq. Q); 

2. solutions with —3 < re < —2 (between the thick curves) interpolate between the Charap e 
Duff classic solutions; 

3. solutions with different values of re do not intersect. 

In what follows, we show that the solutions of interest are indeed those with — 3 < re < 2. In this 
case, the corresponding 4> K does give rise to pairs (4>,a) (cf Eq. ([3])) associated with solutions that 
interpolate between Charap and Duff's instantons with actions I and 2. We also show that item 3 
above is true as long as re < —2, which is precisely the condition for a(r), as defined by Eq. ([6]), to 
be real- valued. 



B. Global behavior of the solutions 



We start by noting a useful local property of the solutions (f) K . Let k±,K2 G K. If k\ < K2, then, 
for r — 2m <C 1: 

^iW<^W. ( 12a ) 

0« 1 (r)<^ 2 (r), (12b) 

which follows directly from Eq. (jlip. Therefore plots of 4> K (r) (as in Fig. [T|) corresponding to 
different values of k do not intersect for r sufficiently close to 2m. A careful analysis of the 
ordinary differential equation (ODE) (J5|) leads to a much stronger statement. In fact, we show in 
the Appendix that for K\ < K2 < —2: 

K1 (r) < cp K2 (r) We(2m,oo), (13a) 
0'« 1 (r)<C(r) Vr€(2m,oo). (13b) 

It follows from Eq. (|13ap that plots of (ft/sll) corresponding to different values of k do not intersect 
for any r G (2m, 00) as long as At < —2 [22|. We also note that taking K2 = —2 in Eq. (|13b|) leads 
to </4( r ) — l/ 7 " 2 whenever k < —2, which is precisely the condition needed to guarantee that a is 
real. 

Consider now a fixed value of k G (—3, —2). A direct consequence of the above discussion is 
that the associated action S K is finite in this case. Indeed, using Eq. (j!3bj) with k\ = —3, we see 
that < 4 > ' K { r )i an d hence <f> K is a strictly increasing function. Moreover, it follows directly from 
Eq. (fTHa]) that 

-%<Mr)<xz-i> for-3< K <-2, 
and thus 4> K (^) is limited. Therefore, the limit 

C K := lim Mr) (14) 

r— »oo 

exists and C K G (0, ^] for —3 < k < —2. On the other hand, it is easy to show, directly from 
Eq. ©, that the asymptotic behavior of <j> K (r), also for —3 < k < —2, is given by 

4>n{r)^C K -- + \ e -^ (r»l) (15) 
r £U K 

for some A G M, which implies a — > 0. This shows that the action can be calculated from Eq. Q 
and that it is indeed finite. It is also clear from Eqs. (jlOap and (|15p that these solutions have unit 
electric charge and therefore correspond to dyons. 

Figure [2] shows a representative subset of the family of solutions (f> K , with k £ [—3,-2] (these 
plots were generated as those in Fig. [T]). Notice that it numerically illustrates all the analytical 
results obtained above. 

It follows from Eq. ([Hj) that, given R > 2m, -gy < (f) K (oo) - (f> K (R) < ^. Thus 

TTl 

MR) + -R2< Moo) < MR) + ^ 

and this leads to accurate estimates for the action via Eq. ([9]): 

S K = 1 + 4m 6(00). 




FIG. 3: Action S K (associated with the solution <f> K ) as a function of k. 
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In fact, data from Fig. [2] yield estimates of S K with an error of one part in 10 6 . We note that the 
above expression, on account of Eq. (|13a|) . limits S K to 1 < S K < 2 if k € [—3, —2]. 

Figure [3] shows how S K depends on the parameter n for —3 < k < —2. We see that the action 
varies continuously with the parameter k. It is worth noting that the derivative of S K becomes 
singular exactly when the physical character of the solution abruptly changes from a dyon to a 
monopole, at k = — 3. 



IV. ANALYTICAL EXPRESSION FOR THE SOLUTION 

So far we have proved the existence of the solutions and presented them numerically. Now 



we make use of the mapping method introduced in [ll|, [12J to obtain an analytical expression for 
the solution defined for all r S [2m, oo). This method is based on the fact that, if the solution 
exists and has no singularities in its domain, then there must be an angular sector of the complex 
plane containing the domain in which the function is analytic. If we compactify this region to 
the complex unit disk, making a transformation from z to a new variable u, then there should be 
a power series expansion for the solution in uj which is convergent for |u;| < 1. To employ this 
method, first we use the shifted variable s = r — 2m. Then we define a new variable related to the 
old one (we are considering the analytic continuation of the solution in the complex plane) by the 
following conformal transformation 11] 



(l + s /fl)V°-l 

S — > W = ; —, , 16) 

(l + s/^) 1/a + l 

where R and a are chosen in order to avoid the first singularity of the analytic continuation of </>(s) 
in the complex plane. We found that the choice R = m and a = 1 apparently gives an infinite 
radius of convergence for the power series (evidence for this will be presented in what follows). The 
inverse transformation is readily obtained 



w — > s = R 



1 i \ a 
1 + to 



— LO 

Using this change of variables in Eq. ([5]) yields 



1 



(17) 



(w - l) 2 ^^ + 4rra/^ - 2^ = 0, (18) 



■ !-—r + Amip — 

du> z cIlo 

where ip(ui) = </>(s(u>)+2m). In this new variable the Charap and Duff solutions become polynomials 
of degree 1 and 2. The family of solutions considered here thus interpolates, in the new variable, 
between a straight line and a parabola. 

We want to write the solutions to Eq. (|18p as a power series in u. As we have already 
discussed, we are interested in solutions with <j)(2m) = ip(0) = Therefore, if we write 

ip{oj) = — Y^=q ^w n+A , we have A = and bo = — |. Substituting the series expression in 
Eq. (|18p we see that b\ = ^, 62 is left free and 63 = 0. The relationship between 62 arid k is given 
by &2 = ^x 2 "- All the other coefficients 6 n +i with n > 3 can be determined via the relation 



2b n [n(n - 1) + 1] - 6 n _!(n - l)(n- 2) -4^(9 + l)6 n _ 9 6 3+1 

b n +l — ]~2 7T 

(n z — 1) 

After evaluating the parameters b n we can go back to the original variable r: 

1 ~ , / 2m 



(19) 



m ^— ' ' \ 

n=0 
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where b\ depends on k via 62- The first few terms of (I20p are displayed below: 



1 



2m 



Mr) = --r- + 7T- 1 + 

4m 2m \ 

k 2 + 5k + 6 / 2m \ 
15m \ r J 



K + 2 



4m 

5 

) +0 



2m V 



r y 

2m V 
r y 



k 2 + 5k + 6 
16?n 



2m\ 4 

-) + 




FIG. 4: Plots of the series expressions obtained before (solid line) and after (dotted line) applying the 
mapping method. The dashed line corresponds to the numerical solution obtained in the previous section. 
Notice that the series solution given by the mapping method (truncated here at n=30 000) is virtually 
indistinguishable from the numerical solution. In all the 



Figure [J] shows a typical solution obtained using this procedure and compares it with its coun- 
terparts obtained by the methods of the previous section. It illustrates the fact that the series 
solution whose first terms are given by Eq. (jlip has a radius of convergence not greater than 2m 
and that the series solution obtained by the mapping method has an apparently infinite radius 
of convergence (see below). It also shows that the latter is virtually indistinguishable from the 
numerical solution obtained in the previous section. 

Figure [5] shows the behavior of b n for a typical value of k. We observe numerically that \b n \ — > 
as n — > 00 so that \b n \ is certainly limited by a constant c > 0. On the other hand, the series 
J^^Lq cuf 1 (whose sum is c/ (1 — to)) is absolutely convergent for |a>| < 1 and using the comparison 
test we see that, as far as the above numerical argument is justified, the series solution will be 
absolutely convergent for cj £ [0, 1), i.e., for any r extending from 2m to infinity 231 ] . 



V. CONCLUDING REMARKS 

We have studied, in a systematic way, smooth spherically symmetric self-dual solutions of Yang- 
Mills theory in the Euclidean Schwarzschild background. We showed that our approach recovers 
all the previously known solutions and leads to a new one-parameter family of instantons with 
continuous actions in the range 1 < S < 2. After studying the global behavior of our solutions, we 



exhibited them numerically and employed the mapping method Hi, LL2|] to express them analytically. 
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FIG. 5: Plot of the coefficients b n as a function of n for n = —2.5. 



Finally, we examined the physical properties of this family and showed that it consists of dyons of 
unit (magnetic and electric) charge which interpolate between the Charap and Duff's monopole and 
dyon. Furthermore, the singularity in the derivative of the action S K with respect to the parameter 
k (Fig. [3]) signals the abrupt change in the physical character of the solution (from dyons to a 
monopole). 

It is interesting to note that it has been argued, on physical grounds, that no time-dependent 
solutions (modulo gauge transformations) exist in Euclidean Schwarzschild space [lo| . As a result, 
our approach would exhaust all possible smooth spherically symmetric solutions (static or not) in 
Schwarzschild space. 

We end with a more mathematical note. Recent results impose strict constraints on the action 



spectrum of instantons defined on ALF geometries 13], of which the Euclidean Schwarzschild space 
is an example. In particular, it is possible to show that an instanton defined on the Euclidean 
Schwarzschild geometry will necessarily have integer action if it satisfies a certain rapidly decaying 
condition, which means that the curvature decays faster than 1/r 2 as r — > oo. We note that the 
noninteger solutions presented in this paper do not satisfy this assumption, see Eq (1151), s o no 
contradiction arises here. On the contrary, this discussion shows that the hypotheses of [13|] are, 
in a sense, as weak as possible. 
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APPENDIX A: GLOBAL ANALYSIS OF THE SOLUTIONS 



We show here that Eqs. (|13[) hold whenever K\ < K2 < —2. We start by noting that, if k < —2, 
then 

M < ^ = <f>L 2 (r) Vr€(2m,oo). (Al) 
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This condition is clearly true for r — 2m <C 1 (see Eq. (|12ap ). Moreover, if this condition ceased to 
hold there would be a ro such that cj)' (ro) = \- But then, by the uniqueness theorem for ODEs, 
we would have 4> K (r) = c — 1/r for some c which, together with the condition on r = 2m, leads to 
= <P-2(f) and therefore k = —2, contradicting our assumption. 
Now we prove that, if k± < K2 < —2, then 

K( r ) < ^(r) Vr€(2?n,oo). (A2) 

As already discussed, this holds for r — 2m <C 1 (see Eq. (|12b|) ). Suppose that there exists 
r G (2m, oo) such that <f>' K1 (r) = <ft' K2 (r) and let ro be the least value of r with this property. Since 
<j)' K is continuous (due to the general theory of ODEs), we have 

4W<i(r) VrG(2m,r ). (A3) 
Integrating both sides of (|A3p from 2m to r G (2m, ro] then leads to 

<M0 < <f> K2 (r) Vr G (2m, r ], (A4) 
since 4> K (2m) = — l/4m for all k. On the other hand, Eq. (JSJ yields 

i(r - 2m)(<4 - <)(r ) + r ^(r ) - ^) (<^ 2 - w )(r o ) = 0, (A5) 

since we assumed that — </>' K2 )(ro) = 0. It follows from Eqs. (|A1[) and (|A4j) that (<^" 2 — )(ro) > 
0, i.e., J^(0K2 ~~ 0Ki)( r o) > 0- I n this way, 4>' K2 {tq — e) < <f>' Kl (ro — e) for e positive and sufficiently 
small. But this contradicts Eq. (|A4|) . Therefore, no such ro exists and Eq. (|A2|) follows. 
Finally, integrating Eq. (IA2|) from 2m to r we see that, for k\ < K2 < —2, 

K1 (r) < </> K2 (r) Vr G (2m, oo). (A6) 

This concludes the proof. 
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